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Abstract 

We consider a soluble model of large c/> 4 -graphs randomly embedded in one compactified dimension; 
namely the large-order behaviour of finite-temperature perturbation theory for the partition function 
of the anharmonic oscillator. We solve the model using semi-classical methods and demonstrate the 
existence of a critical temperature at which the system undergoes a second-order phase transition 
from D = 1 to D = behaviour. Non-trivial windings of the closed loops in a graph around the 
compactified time direction are interpreted as vortices. The critical point has a natural interpretation 
as the temperature at which these vortices condense and disorder the system. We show that the 
vortex density increases rapidly in the critical region indicating the breakdown of the dilute vortex gas 
approximation at this point. We discuss the relation of this phenomenon to the Berezinskii-Kosterlitz- 
Thouless transition in the D = 1 matrix model formulated on a circle. 



1 Introduction 



Models of random graphs have applications to theories of polymers and random surfaces. An 
example of particular interest is the D — 1 Hermitian matrix model pi ^ |3|, [3j formulated on 
a circle of radius R which is equivalent to matrix quantum mechanics at a finite temperature 
ksT = 1/27T.R. The perturbative expansion of the matrix path-integral for the case of a quartic 
potential generates a sum over 4 Feynman graphs, each one dual to the quadrangulation of a 
random surface; these surfaces being randomly embedded in S 1 . In the double-scaling limit, the 
resulting continuum theory can be interpreted as an XY model coupled to two-dimensional quan- 
tum gravity or as non-critical string theory in one compactified embedding dimension (an extra 
dimension corresponding to the Liouville mode also arises). An important difference between a 
model of graphs embedded on the circle and one formulated on the real line, is that the former ad- 
mits vortex configurations where a closed loop in the graph wraps around the compactified time 
direction. Gross and Klebanov have argued that these vortices are suppressed at large radius 
but, at critical value of the radius R c , the system undergoes a Berezinskii-Kosterlitz-Thouless 
(BKT) U transition to a high-temperature phase in which the vortices condense and disorder the 
embedding coordinate. They find that the embedding coordinate is completely randomized for 
R < R c and effectively decouples leading to zero-dimensional behaviour in the high-temperature 
phase. 

Unfortunately the D = 1 matrix model on a circle cannot be solved in the double scaling limit 
by the usual methods which are effective in the case of a non-compact target space. The reason 
is that the integral over the angular degrees of freedom of the matrix field, which can be done 
exactly on the real line to yield a quantum mechanical problem involving the eigenvalues only, is 
no longer tractable on the circle. Indeed this difficulty seems to be precisely due to the presence 
of vortices in the latter model: at low temperature the model can be solved approximately by 
considering only the U (N) singlet sector, however this appears to be equivalent to neglecting the 
vortex contribution FL 

Similar vortex configurations occur in any model of random graphs embedded on a circle. In 
general, these graphs need not correspond directly to discretized random surfaces as they do for 
the matrix model. In this letter we present such a model which exhibits vortex condensation 
accompanied by an abrupt transition from D = 1 to D = behaviour and has the additional 
virtue of being exactly soluble. The model involves an ensemble of large 4 -graphs randomly 
embedded in a single compactified dimension and is equivalent to the large-order behaviour of 
finite-temperature perturbation theory for the anharmonic oscillator. The graphs we consider 
carry no index structure, which means that there is no unique prescription for attaching plac- 
quettes and thereby defining a discretized surface. Thus, it will not be possible to draw any 
definite conclusions for one-dimensional string theory from our results. Nevertheless, the model 
exhibits a non-trivial critical behaviour which is similar to that expected at the BKT transition 
point of the D — 1 matrix model and therefore merits study. 

1 Boulatov and Kazakov Q have considered the adjoint sector of the theory and have shown that this corre- 
sponds to the contribution of a single vortex-antivortex pair. 
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The large-order behaviour of perturbation theory for the quantum anharmonic oscillator de- 
fined by the Hamiltonian, 

H=(p 2 + ^)/2 + g<p\ (1.1) 

is determined by the instanton solution |7|, || , which is the classical motion of a particle in the 
inverted potential V = —x 2 /2 + x 4 /4 (see Figure 1). At finite temperature T, the relevant 
classical motion is constrained to have period f3 = 1/ksT. It is clear that there is a minimum 
period for which there exists a non-trivial instanton solution x(t). The period is that of an 
infinitessimal simple harmonic motion about the quadratic minimum of V and is determined 
by the curvature of the potential well at that point. Below this critical period, the large-order 
behaviour of perturbation theory is governed by the trivial solutions x(t) = ±1. Although this 
phenomenon has no special signifigance for the anharmonic oscillator itself 0, we will argue that 
it has a natural interpretation as a second-order phase transition in the corresponding model of 
randomly embedded graphs. In particular, we will show that the specific heat is sharply peaked 
at the critical temperature and exhibits a finite discontinuity at this point. In addition we will 
show that density of free vortices undergoes a crossover in the vicinity of the critical point from 
a low temperature regime of tightly bound vortex-antivortex pairs to a high temperature regime 
populated by free (anti-)vortices. Thus we will interpret the critical temperature as the point at 
which vortices condense and disorder the system. 

In Section 2, we define the partition function of the model as a sum over Feynman graphs with 
a fixed large number of vertices and identify the vortex configurations. In Section 3 we solve the 
model using semi-classical methods as described above and exhibit the critical behaviour. In the 
process, we give some new analytic results for the large-order behaviour of perturbation theory 
for the anharmonic oscillator at finite temperature. Section 4 is devoted to a discussion of these 
results. 



2 The model 

The partition function for the anharmonic oscillator defined by the Hamiltonian ( |1 . 1|) , at a 
finite temperature T, is given by 



Z((3,g) = Tr 





= J £>0exp 


["/„ 









(2.2) 



with f3 = l/ksT. Finite-temperature perturbation theory is generated by expanding the partition 
function as an asymptotic series in the coupling constant g, 



Z((3,g)~J2Z k ((3)g h 

k=0 



(2.3) 



The coefficients are given by the sum of all connected <j) 4 vacuum graphs G with k vertices 
and the Feynman rules instruct us to assign a time coordinate £, £ [0,/3] to each vertex and a 

2 The occurence of the minimum period was noted in this context in H. 
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finite-temperature propagator, 



CO 



(2.4) 



m=— oo 



to each link (ij) and then integrate over each tj. The propagator has a geometrical inter- 
pretation as a sum over the topologically inequivalent ways of embedding the link (ij) in the S l 
obtained by identifying the endpoints of the time interval [0,/3] (see Figure 2). Each embedding 
of a link (ij), is weighted by a factor exp(— i/3) in the summation ( |2.4|) , where / is the length (in 
units of /3) of the image of (ij) in the target space. 

Applying the Feynman rules and summing over all graphs with k vertices we have, 



For each graph G, a choice of the time coordinate {tj} for each vertex and of a winding number 
rriij for each link specifies an embedding of G in S . The coefficient Zk is a sum/integral over all 
such embeddings of all /c-vertex vacuum graphs G, each one weighted by the symmetry factor 
S(G) and by exp (-L/3) where L is the total length of the image of G in the target space. In 
what follows, we will consider Zk(/3) itself, for fixed large k, as defining the partition function 
of a theory of graphs randomly embedded in S 1 , where the "action" of a "configuration" (ie a 
specific embedding) is just its target-space length. From this point of view, k — > oo is the bulk 
limit of the model, and we can define the usual thermodynamic quantities. In this limit, the 
number of graphs grows as ~ (16) fe (A; — 1)! and we define a normalized partition function 
Z^ = Zk/rik which can be thought of as the contribution of a typical graph. The free energy per 
vertex is defined as, 



and the corresponding specific heat is given by c v = —f3 2 d 2 f/d(3 2 . In Section 3 we will give the 
exact solution of the model in the thermodynamic limit by applying the standard semi-classical 
analysis of large orders in perturbation theory to this finite temperature case. In particular, we 
will show that / and c v remain finite as k — > oo and that the latter exhibits a discontinuity as 
a function of temperature indicating a second-order phase transition. It is convenient to express 
the partition function, Zk((3), as a Laplace transform, 



The inverse Laplace transform, Zk(L), is the corresponding microcanonical partition function for 
configurations of fixed target-space length L. 

For each configuration, it is possible to assign a vortex number me to every closed loop I in 
G. Specifically, each term obtained by expanding the product of propagators, 




(2.5) 



/ = - lim -log|Z fc (^)| 



(2.6) 




(2.7) 



CO 



J| exp[-\ti-tj + mij/3\] 




(ij)&£rriij=-QO 
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is characterized by a vortex number me = J2(ij)£e m ij which labels the homotopy class of the 
corresponding embedding of I in S 1 . Negative values of mi correspond to anti- vortices and, by 
analogy with the two-dimensional XY-model, we interpret loops for which me = as tightly- 
bound vortex-antivortex pairs. The total vortex number of a configuration is given by summing 
over all closed loops in G, M = J2e m e- Because the sum over configurations ( |2.5| ) is symmetric 
under m^- — > —mij, vortices and anti- vortices occur with equal statistical weight and therefore 
(M) = 0. 

For low temperatures, the summation for each propagator in ( |2.5p is dominated by the m^ = 
term. An isolated (anti-)vortex necessitates at least one of the my being non-zero and is thus 
suppressed by a factor of exp (— (3). Thus, at low temperature, we expect that the links with 
non-zero m^ in a typical configuration are well seperated and therefore that most vortices and 
antivortices are tightly bound in pairs. The unbinding of a vortex-antivortex pair is accompanied 
by an increase in the target-space length of the configuration by two units. Thus, at least in 
this dilute low-temperature ground state, the expected number of free (anti-) vortices increases 
linearly with the expectation value of the quantity L. In addition, when the corresponding length 
per vertex p = (L)/k becomes 0(1), it follows that the free vortices have become dense and the 
dilute gas approximation is no longer valid. In this sense we will refer to p as a "vortex density" . 
For any given configuration, L is also the corresponding action and so the vortex density in the 
thermodynamic limit is given by the internal energy per vertex, p = df / df3 (this formula can be 
obtained by differentiating ( |2.7| )). In the next section, we calculate p explicitly and show that 
the system leaves the dilute regime in the vicinity of the critical temperature. This suggests that 
the phase transition corresponds to a disordering of the system by the free vortices. 

3 Semi-classical solution for k —> oo 

In this section we follow the standard approach to the large order behaviour of perturbation 
theory for the anharmonic oscillator |7| (for a pedagogical exposition see P|). Although this 
problem is frequently formulated at finite temperature as a way of regulating the infra-red diver- 
gences which occur at zero temperature, the explicit evaluation of the temperature dependence of 
large-order behaviour given below appears to be new. The coefficients Zk(f5) in ( |2.3| ) can be ex- 
pressed in terms of the imaginary part of the original partition function (|2.2|) for the anharmonic 
oscillator as, 

reflecting the existence of a cut in the complex (/-plane along the negative real axis. For k — > oo 
the integral is dominated by the g — > region of the integrand. In this region, the functional 
integral (|2.2| ) for Z(g) can be evaluated accurately in the saddle-point approximation. The saddle- 
point solutions, x{t) = 2^—g(j)(t), are given by the appropriate periodic classical trajectories of 
a particle moving in the potential, V = —x 2 /2 + a; 4 /4 with turning points x_ < x + and total 
energy E = V(x±) (see Figure 1). At finite temperature the relevant trajectories have period (3. 
As discussed in the Introduction, there is a minimum period (3 C for which a non-trivial solution 
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exists. Expanding the potential in the vicinity of the positive minimum x c = 1 we find, 

V = ~ + X -u 2 (x - x c f + 0((x - x c ) 3 ) (3.10) 

with uj = y/2. Thus as the total energy of the particle approaches a critical value E c = —1/4, 
its motion becomes an infinitessimal harmonic oscillation about the well-bottom with period 
C = 2ix/uj = \^2ti. Although n-instanton corrections corresponding to paths with periods 0/n 
can occur at low temperature, these solutions do not occur for < 2/3 c . 

The relevant instanton solutions, x(t) satisfy the classical equation of motion, 

x- x + x 3 = (3.11) 

with boundary conditions x(0) = x(0). For > C , there are two one-parameter families of 
such solutions generated by time-translations and reflections, x — *> —x. For definiteness, we will 
consider the particular solution such that x(0) > is the turning point of the motion x_. This 
solution is given as a Jacobian elliptic function §|, 

x(t) = x_dn k (3.12) 

with imaginary modulus, k = \Jl — x\jx 2 _. 

The resulting saddle-point expression for ImZ(g) depends on (3 as 0, 

where 1(0) is the classical action of the trajectory, which is conveniently determined by the 
relation 

together with the condition I (0c) — /?c/16 for the action of the limiting infinitessimal harmonic 
motion. The pre-exponential factor in (|3.13|) is given by the determinant of gaussian fluctuations 
about the saddle-point and is proportional to 1/ ' ^fg reflecting the contribution of the zero- mode 
corresponding to time-translation of the instanton solution. Adopting the convention of ||, all 
functional determinants considered (see also ( 3.22| ) below) are normalized with respect to the 
zero coupling result. 



The formulae ( |3.14| , |3TT3"| ) demand that we find the dependence of the turning point energy 
E = V(x±) on the period 0. In particular we are interested in the behaviour of E(0) and its 
derivatives in the critical region (0 — C ) << 1. The inverse function, 0(E), can be expressed as 
a first integral of the classical equation of motion ( 3.11]) , 



„ = 2 i £ = t E dv w^^ (3 . 15 ) 

<2{E - V(x)) J e c dV / 2 (E - V) 
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where W(E) = x + — X- is the width of the well at energy E < 0. The above relation (Abel's 



integral equation) can be solved for W by noting that the second integral in Q3.15 ) is propor- 
tional to the Riemann-Liouville integral^ [|I(J which defines the fractional derivative D~z. In 
particular, using the composition rule for fractional derivatives [Kj D 6 D 1 = _D 5+7 , we have 
(3{E) oc (D~2DW)(E) = {D^W){E), which implies that W(E) oc (zH/3)(£) or, reintroducing 
numerical factors, 

W(E) = - [ E dE' . (3.16) 
7T Je c . o(F, - FA 



*2(E - E> 

However, in the present case of a quartic potential, W(E) is given explicitly by, 



W(E) = yl + yjE — E c — yl — \JE — E c 

- (3-17) 
Comparing ( |3.16| ) and ( |3.17| ) we find that (3{E) is analytic for all E < and has a Taylor series, 



P(E) = & + p c £ a n (E - E c ) n (3.18) 



with coefficients, 



71=1 



Standard theorems on implicit functions imply that E is an infinitely differentiable function of 
(3 in some neighbourhood of f3 c . The Taylor series for E{(3) in powers of {(3 — /3 C ) can be derived 
by inverting the series ( [3.18 ). Although the general term in this series cannot be expressed in 



closed form, the coefficients can be generated systematically using the algorithm described in the 
appendix. For the present purpose, we will require only the leading terms, 

EW) = -1 + ±.{p - fi c ) + 0(((3 - &) 2 ) (3.20) 

For (3 —* oo, E{(3) has a double expansion in powers of (3 and e~^. A straightforward calculation 
yields, 

E{(3) = -l6e~P - 64(3/5 - 10)e" 2/3 + 0(/?V 3/3 ) (3.21) 
The first term in this series yields the standard zero-temperature result || when substituted 



for E in (|3.13|). We computed E{(3) numerically for the whole range (3 C < (3 < oo by solving 



equation (|3.15| ), the resulting graph is shown in Figure 3. We also computed I {(3) by integrating 



our numerical solution for E((3). 

3 This integral can also be thought of as arising from the standard Duhamel integral jll| by explicit 
differentiation. 
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For (3 < /3 C , the only solutions of the equation of motion are the trivial ones x(t) = ±1. The 
corresponding contribution to the partition function is, 



imz( ^KSf) exp 



p 



16a 



(3.22) 



where the determinant prefactor is independent of g as the trivial saddle-point has no zero- 
mode. Applying equation ( |3.9|) , we find that (3 C separates two "phases" characterized by different 
asymptotic behaviour of Zk- For (3 > /3 C we have, 



; 2 \2tt 8(3 ) 



(3.23) 



whereas for (3 < (3 C 



:-i) fc 2 



V2sinh/3 
, sinh y/20 




F(k) 



(3.24) 



In the low-temperature phase, Zk exhibits the same I~ k k l l 2 {k— 1)! growth, as the zero-temperature 
result. In contrast, the large-order behaviour of perturbation theory in the high-temperature 
phase is identical (up to a /c-independent prefactor) to that of the perturbative expansion of the 
ordinary integral, 

P (y 2 



dyexp 



4 



+ gy H 



(3.25) 



In this sense, the large-order behaviour of perturbation theory exhibits a transition from D = 1 
to D = behaviour at (3 = (3 C . 



Equations ( |3.23|) and ( 3.24 ) constitute an exact solution for the partition function of the model 
of random graphs considered in Section 2 in the thermodynamic limit. Although the partition 
function itself grows wildly with k reflecting the k\ growth in the number of Feynman diagrams, 
the free energy per vertex of a typical graph fl2.6|) remains finite as does the corresponding specific 
heat. For (3 > /3 c we find, 
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BE 



1 + E 2 ((3) 



8(3 1{P) 4P((3) 



(3.26) 



while for f3 < (3 C , c v = 1. Using equations ( |3.14| , |3.20| , |3.26|) we find that c v — > 19/3 as (3 { (3 C . 
The resulting discontinuity in the specific heat at (3 = (3 C implies that the system undergoes a 
second-order phase transition at this point. The result of a numerical evaluation of c v is shown 
in Figure 4. As expected for a second-order phase transition, the specific heat has a narrow peak 
in the critical region. 
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In Section 2 we argued that the density of free vortices is given by the expectation value of 
the target-space length per vertex L/k which is exactly the internal energy df /d(3. Hence, for 
(3 > j3 c , we have p = —E((3)/AI((3) while for f3 < f3 c , p — 1/(3. The corresponding graph of the 
vortex density against temperature is shown in Figure 5. At low temperature, the vortex density 
is exponentially suppressed, 

p ~ exp (3.27) 



This implies that the low-temperature ground state of the system can be thought of as a dilute 
gas of tightly-bound vortex-antivortex pairs. In the high-temperature phase the vortex density 
is 0(1) and increases linearly with temperature indicating a breakdown of the dilute gas approx- 
imation. The numerical results shown in Figure 5 reveal a rapid crossover between these two 
regimes in a narrow region near the critical temperature k B T c = l/v^2vr. In particular, the vortex 
density increases by a factor of about 6.35 in the interval 0.75T C < T < T c . This supports an 
interpretation of the critical temperature as the point at which vortex-antivortex pairs unbind 
and disorder the system. There is a qualitative similarity between our numerical results for p 
and an evaluation of the vortex density in a Monte-Carlo simulation of the two-dimensional XY 



model presented in [12 



4 Conclusions 

In this paper we have presented a novel interpretation for the large-order behaviour of finite- 
temperature perturbation theory. In particular we have shown that a fixed large order in the 
perturbation series for the partition function of an anharmonic oscillator defines a non-trivial 
statistical mechanical model with a sensible thermodynamic limit. The model exhibits two phases 
separated by a second-order phase transition and it is important to establish the universality class 
to which this critical behaviour belongs. Although this transition has the same interpretation as 
the BKT transition in the two-dimensional XY model, the latter is a transition of infinite order 
having no discontinuities in the derivatives of the free energy. 

The behaviour of Z^ for (3 < (3 C is identical to the large-order behaviour of the perturbative 
expansion of an ordinary integral. The perturbation series for the integral Z* is just a sum over 
all 4 graphs with no embedding dimension. This suggests that, in the high-temperature phase 
of the model, the embedding coordinate is completely randomized by the gas of free vortices and 
decouples leaving an effectively zero-dimensional theory. This is exactly the behaviour expected 
at the BKT transition point of the D — 1 matrix model p[. However, due to the unconventional 
procedure of summing over all Feynman graphs with equal weight irrespective of genus, there is 
no direct correspondence between our model and the double-scaling limit of the D = 1 matrix 
model. In addition, we note that the estimates for the critical temperature of the matrix model 
given in 0, Q differ from the critical temperature found here by a factor of two. On the other 
hand, the occurence of vortices is a purely local phenomenon on any graph or discretized surface 
and it is expected that the BKT transition should occur independently at each order in the genus 
expansion. Hence one might expect the same critical behaviour to appear even in the unweighted 
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sum over different genera considered here. Thus it is possible that our model exhibits the same 
universal behaviour as the XY model coupled to two-dimensional quantum gravity. 

The model presented in this paper has several unusual features. As in the matrix models, the 
sum over configurations is generated by the perturbative expansion of a functional integral and 
can be evaluted using semi-classical methods. However, in our case, the thermodynamic limit of 
large graphs is taken explicitly rather than by tuning the coupling constant to a critical value 
(indeed the latter procedure is not available here as the perturbation series has zero radius of 
convergence). In principle we can obtain a continuum theory only at the critical temperature 
/3 = f3 c . Finally we note that the model considered here can be generalised to any theory where 
the large-order behaviour of perturbation theory may be evaluated by saddle-point methods. 

The authors acknowledge useful discussions with Jan Ambjorn, Paul Ginsparg, Tomek Jaroszewicz, 
and Ian Kogan. PSK thanks the Aspen Center for Physics, where part of this work was com- 
pleted, for its hospitality. 
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Appendix 

In this appendix we give a general formula for the inversion of a power series of the form, 

oo 

w = /(*) = £ a nz n (4.28) 

n=0 



which provides an algorithm for generating the coefficients of the series (|3.20|) for E(f3). We 
consider the inversion of ( |4.28| ) in an open neighbourhood of z = 0; accordingly we demand that 
f'(0) 7^ and set /(0) = Wq. The inverse series for z(w) can be written as, 



Z=-(W-W )+Y, ^ „2m-l -£•■■ £ ^k-m+1,0) C h . . . C lm (4.29) 

1 h=0 l m =0 



m=2 



nil 



a 



where 



a 



Cl2 Oi 

O3 &2 CL± 

di a^i a/_2 

ai+i ai a/_i 






a 1 
a 2 



(4.30) 



More explicitly, if 



w = Wq + az + bz 2 + cz 3 + dz 4 + . . . 



(4.31) 



flPgp gives 
1 , 



: -( w-wo)-—(w-wo) 2 -^(ac-2b 2 )(w-wo) 3 -\(a 2 d~5abc + 5b 3 )(w-w )' i + . ■ ■ (4.32) 
a a 6 a 5 a 7 



Figure Captions 

Figure 1. The potential V(x). 



Figure 2. The geometrical interpretation of the finite-T propagator (|2.4Q as a sum over topo- 
logically inequivalent embeddings of the link (ij) in S 1 . 



Figure 3. A graph of E as a function of (3/ (3 C 
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Figure 4. A graph of c v as a function of k B T = 1/(3. 
Figure 5. A graph of p as a function of ksT. 
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